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ABSTRACT 

The Gondwanian-derived Asian pomatiopsid radiation is taxonomically complex, diversity-rich, 
and widely deployed geographically. This Asian branch of the family has coevolved with such 
human trematode parasites as Schistosoma and Paragonimus: it is ideally suitable for studying 
patterns and processes of evolution over 1 00 million years. Cytochrome c oxidase subunit I gene 
sequences are used here to elucidate taxonomic relationships from the subspecies to familial 
level. In Chinese literature, pomatiopsid taxa have been classified in the Hydrobiidae; what are 
the genetic relationships between Hydrobia and allied taxa classiried as pomatiopsid? 

Sixteen sequences, ranging in length from 578 to 645 nucleotides, are aligned from 11 species 
of nine genera assigned to seven families, four of which are rissoacean. Five different phyloge- 
netic analyses are concordant: (1) the pomatiopsid taxa are in one distinct clade, the other ris- 
soaceans form a second clade; (2) truncatellids are more closely allied to the hydrobiids than to 
the pomatiopsids; (3) the rissoid Setia is part of the truncatellid-hydrobiid clade; (4) two sub- 
species of Oncomelania are clearly divergent; (5) triculine taxa appear divergent from poma- 
tiopsine taxa. However, the Tricula sp. node is weakly supported. 

Individuals of a population differ by an average of 0.005 n 0.004 nucleotide differences/site; 
the subspecies of Oncomelania differ by 0.148 0.004; the two species of Hydrobia differ by 

0.162 (range of 0.161 - 0.163); the triculine genera Tricula and Gammatricula differ by 0.132 
(range of 0.130 - 0.133); the pomatiopsid subfamilies Pomatiopsinae and Triculinae differ by 
0.179 ± 0.020; the families Hydrobiidae and Pomatiopsidae differ by 0.267 ± 0.016. Non-ris- 
soacean and rissoacean taxa differ by 0.274 ± 0.023. 

Key words: systematics, cytochrome c oxidase subunit I, COI, gene sequences, phylogeny, 
Rissoacea, Pomatiopsidae, Hydrobiidae, Truncatellidae, Rissoidae, China. Jamaica, Bulgaria, 
Denmark. 


INTRODUCTION 

This study is one of a series (reviewed in 
Davis, 1992) aimed at understanding the ori- 
gin and evolution of the freshwater snail fam- 
ily Pomatiopsidae in Asia. There are com- 
pelling reasons to pursue such studies: (1) 
The family is ideally suited for in-depth studies 
of biogeography and evolution. The family is 
of Gondwanian origin with genera found in 
South Africa, South America, northern India, 
and Australia. Davis (1979) hypothesized an 
introduction of early pomatiopsids from the 
northeastern Indian Plate into northern Burma 
and western China with subsequent distribu- 
tion throughout southern China, Japan, and 
the Philippines reaching North America via 


Bering Strait. (2) Pomatiopsids are ideal for 
studying patterns and processes of evolution 
over the past 100 million years. The family is 
taxon rich and widely deployed geographi- 
cally. A series of anatomical studies on poma- 
tiopsid taxa have yielded a rich database of 
characters and character-states enabling the 
establishment of testable phylogenetic hy- 
potheses of the evolution of the family. Each 
phylogeny is tested with the addition of data 
from a newer study. The latest phylogeny 
(Davis, 1992) has not falsified the previous 
phylogenies. These phylogenies are mapped 
on area cladograms testing and reinforcing 
the biogeographic hypothesis. (3) A diverse 
array of Asian pomatiopsids are important in- 
termediate hosts of human trematode para- 
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sites. These studies are essential to study co- 
evolution of the blood parasite Schistosoma 
with prosobranch snails. Three genera of the 
Pomatiopsidae — Oncomelania, Neotricula 
and Robertsiella — are involved in the trans- 
mission of the human blood fluke Schisto- 
soma In China, the lower Mekong River, and 
Malaysia. The schistosome occurring in 
China is also distributed in Japan, Taiwan, the 
Philippines and Sulawesi. The genus Schisto- 
soma is, as are the Pomatiopsidae, of Gond- 
wanian origin; one branch of the Schistosoma 
clade has coevolved with Asian Pomatiopsi- 
dae, with genera distributed in two subfami- 
lies, the Pomatiopsinae and Triculinae (re- 
viewed In Davis, 1980, 1992). (4) Molecular 
data are essential to test the phylogenetic re- 
sults based on anatomy. Allozyme data have, 
to the limited degree they have been applied, 
reinforced the phylogenetic hypotheses on 
the relationships of key genera in the afore- 
mentioned subfamilies (Davis et al., 1994); 
that is, the phylogenies based on allozymes 
and anatomical data are congruent. They 
have been useful to clarify subspecific status 
of populations of Oncomelania hupensis in 
China (Davis et al., 1995). But with the need 
to include numerous taxa for phylogenetic 
analysis, the ability to use allozyme data ef- 
fectively decreases. We have found that mito- 
chondrial gene sequences are ideally suited 
for testing phylogenetic relationships based 
on anatomical data. We have established that 
cytochrome b sequences are useful to deter- 
mine patterns of divergence at the population 
and subspecies level within the genus On- 
comelania (Spolsky et al., 1996). 

There are a number of questions that this 
study was designed to answer: (1 ) Are the Hy- 
drobiidae and Pomatiopsidae truly divergent 
separate families? From Davis (1979) onward 
we have argued that, on the basis of anatomy 
and patterns of development, the Hydrobiidae 
and Pomatiopsidae are highly divergent. This 
is an important question because up to the 
present, Chinese workers and others around 
the world have insisted on including the Po- 
matiopsidae within the Hydrobiidae and thus 
believe that some hydrobiids transmit schisto- 
somes (Malek & Little, 1971; Brandt, 1974; 
Liu et al., 1974; Liu, 1979; Brown, 1980 (re- 
vised to use Pomatiopsidae in 1994); Kang, 
1984, 1986; Malek, 1985). Davis (1979, 1980, 
1 992) has shown that the Hydrobiidae are not 
found in India, China or southeast Asia; fur- 
ther, no hydrobiid transmits Schistosoma. 

Understanding the patterns of divergence 
of these two families is important to under- 


standing the origin and coevolution of ris- 
soacean snails with Schistosoma. One pur- 
pose of this paper is to present additional ev- 
idence that the two families are distinct and 
divergent. 

(2) On the basis of anatomical data, there 
are two distinct subfamilies of the Pomatiopsi- 
dae; the Pomatiopsinae and Triculinae. Al- 
lozyme data reinforced the confamilial status 
of the two generic groupings but did not un- 
equivocally serve to demonstrate two distinct 
subfamilies as anatomical data did (Davis et 
al., 1994). Would the cytochrome c oxidase 
subunit I (CO!) gene sequences serve to 
clearly demonstrate family and subfamily- 
level generic groupings? 

(3) Would the Truncatellidae (represented 
here by only one species) be more closely re- 
lated to the Hydrobiidae or to the Pomatiopsi- 
dae? Davis (1 979) and Ponder (1988) consid- 
ered the Truncatellidae to be closely related to 
the Pomatiopsidae on the basis of anatomical 
data. Analyses of 28S rRNA sequences in- 
volving one species of Hydrobiidae, nine 
species of two genera of Truncatellidae, and 
one species of Pomatiopsidae (Rosenberg et 
al., 1997) placed the Hydrobiidae as an out- 
group to a cluster consisting of two branches: 
one solitary branch included some of the trun- 
catellid taxa; the second branch subdivided 
into two groups, one including the remaining 
truncatellid taxa, the other the pomatiopsid 
species. What would the COI sequence data 
tell us? 

(4) Would COI data support the conclusion 
based on allozyme data (Davis et al., 1995) 
that Oncomelania hupensis hupensis and O. 
hupensis robertsoni are distinct subspecies? 


METHODS 

Taxa Studied 

Pomatiopsidae: Triculinae: Gammatricula 
chinensis Davis, Liu & Chen, 1 990; Tricula sp. 

Pomatiopsidae: Pomatiopsinae: Oncome- 
lania hupensis hupensis (Gredler, 1881); On- 
comelania hupensis robertsoni Bartsch, 
1946. 

Hydrobiidae: Hydrobiinae: Hydrobiaci. pon- 
tieuxini Radoman, 1973; Hydrobia neglecta 
Muus, 1963. 

Rissoidae: Setia turriculata Monterosato, 
1884. 

Truncatellidae: Truncatella pulchella (Pfeif- 
fer, 1839). 


CYTOCHROME OXIDASE-I-BASED RISSOACEAN PHYLOGENY 


253 


TABLE 1. Localities and collecting information for the specimens studied. 


Taxa: Preparation # 

Localities 

Catalog# 

Pomatiopsidae:Triculinae 
Gammatricula chinensis 
414/415/416 

Tricula sp. 

453/454 

China, Zhejiang Province, Kaiwa Co, 

Tong Cun Town, Bai Keng Village; 
118°15'47"E, 29"00'05"N 
China, Sichuan Province, Dayi County; 

Tian Gong Mia Township; Huang Ba Village; 
117 23'16"E, 30"35'26"N 

ANSP 400351 
ZAMIPM0136 

ANSP 400352 

Pomatiopsidae: Pomatiopsinae 
Oncomelania hupensis 
hupensis 93/96 
Oncomelania hupensis 
robertsoni 45/48 

China, Hubei Province, Han Yang County; 

114°0r01"E, 30"34'08"N 
China, Yunnan Province; Dali City, 

Da Jin Ping, Zi Ran Village; 
100"12W'E, 25'=27'06"N 

ANSP 375731 
CIPD 0349 

Hydrobiidae 
Hydrobia cf. pontieuxini 
346/347/351 
Hydrobia neglecta 
435/436/439 

Bulgaria, 1 km W of Nessebar; 

27"71'73"E, 42 65'99"N 
Denmark, Funen Island, Odense Fjord; 
10 32'E, 55"30'N 

ANSP 400353 
ANSP 400354 

Rissoidae 
Setia turriculata 
474/476/477 

Bulgaria, 1 km W of Nessebar; 
27 7173"E, 42 65'99"N 

ANSP 400355 

Truncatellidae 

Truncatella pulchella 
479-480 

Jamaica, W of Falmouth; 
77"39'46"W, 18 29'46"N 

ANSP 400356 


Four molluscan taxa were used as out- 
groups: (1) the polyplacophoran Katharina tu- 
nicata (Wood, 1815), the sequence for which 
was obtained from GenBank (accession num- 
ber U09810); (2) the cerithiacean Cerithium 
atratum (Born, 1778) (Harasewych et aL, 
1997); (3) the muricacean Stramonita haema- 
stoma (Linnaeus, 1767) (GenBank accession 
number U86330, under the name Thais); and 
(4) the rissoacean Sofia turriculata studied in 
this paper. 

Locality data are given in Table 1 . Oncome- 
lania hupensis hupensis, O. hupensis robert- 
soni and T pulchella were brought to the USA 
alive, G. chinensis and Tricula sp. were pre- 
served in 100% methanol, H. cf. pontieuxini 
and H. neglecta were preserved in 70% 
ethanol, and S. turriculata was frozen. Imme- 
diately prior to isolation of DNA, the living 
specimens of Oncomelania and Truncatella 
were quick-frozen at -80 C. 

DNA Preparation 

The methods used for preparing DNA from 
individual snails were described by Spolsky et 
al. (1996), with the following modifications. Al- 
cohol-preserved specimens of Gammatricula, 


Tricula and Hydrobia were soaked 5 min each 
in two changes of 300 pi of ice-cold exchange 
buffer before being placed in lysis buffer. 
Ethanol precipitation and washing were re- 
peated and the final DNA pellet redissolved in 
25 pi of water. 

The quality of DNA was determined by 
electrophoresis through a 1% agarose gel in 
TBE. DNA concentration was determined 
using a HoeferTKOlOO fluorometer. 

DNA Amplification 

PCR was used to amplify a fragment of the 
mitochondrial COI gene using the primer pair 
COF14 (forward: 5' GGTCAACAAATCATA- 
AAGATATTGG 3') and COR722 (reverse: 5' 
TAAACTTCAGGGTGACCAAAAAAYCA 3'). 
COF14 is identical to primer LCO1490 as 
described by Folmer et al. (1994), while 
COR722 is a modification of Folmer et al. 
(1994) primer HC02198. 

Each PCR reaction mixture, in a total vol- 
ume of 50 pi, contained 20-100 ng of genomic 
DNA, 2.5 units of cloned Pfu DNA polymerase 
(Stratagene), 200 pM of each dNTP, 20-40 
pmol of each primer, 20 pg of BSA, and 5 pi of 
10X Pfu reaction buffer. PCR amplifications 


254 


DAVIS ETAL. 


were performed using an M-J Research PTC- 
100 thermal controller with the following cy- 
cling conditions: initial 1 min 30 sec at 95°C, 
followed by 40 cycles of 1 min at 95°, 1 min 20 
sec at 47°C, and 1 min 1 0 sec incremented by 
1 sec per cycle at 73°C. After a final 5 min at 
74°C, reactions were held at 4°C. 

Amplified DNA products were separated by 
electrophoresis through a 1% low melting 
point agarose gel In TAE buffer. The band cor- 
responding to a fragment of the correct size 
was cut out, and the DNA purified using either 
Mlcrocon-1 00 microconcentrators (Amicon) 
after digesting the agarose with agarase or di- 
rectly with Wizard PCR preps (Promega). 
DNA concentration was determined on a Hoe- 
fer TK01 00 fluorometer. 

Sequencing 

Sequences of the COI fragment were de- 
termined by manual cycle sequencing, using 
the Thermo Sequenase radiolabeled termina- 
tor cycle sequencing kit (Amersham) accord- 
ing to their protocol. 

Each reaction mixture contained 40-60 ng 
of purified PCR product, 4 units of Thermo Se- 
quenase DNA polymerase, 2 pi of reaction 
buffer, 4 pmol of either COF14 or COR722 
primer in a total reaction volume of 20 pL One 
fourth of the reaction mixture was aliquoted to 
each of four dideoxynucleotide-specific termi- 
nation mixes containing the four dNTPs at 7.5 
pM each plus, for each termination, 0.1 5 pmol 
of one of the four ^^P-labelled ddNTPs at a 
specific activity of 1500 C/mmole. Cycling 
conditions consisted of 30 cycles of 60 sec at 
95°C, 60 sec at 51 C and 75 sec at 72°C for 
the forward primer and 30 cycles of 60 sec at 
95 C, 60 sec at 56°C and 75 sec at 72 C for 
the reverse primer. Four pi of stop solution 
were added after cycling and samples were 
stored frozen until ready to use. 

Before loading, the samples were heated 5 
minutes at 75"C and 3 pi of each loaded im- 
mediately on the gel. Reaction products for 
each primer were run on three separate gels: 
2.5 h on an 8% Long Ranger gel (FMC), 4 h 
on a 6% gel, and 7.5 h on a 5% gel at a con- 
stant power of 37.5 Watts, providing complete 
sequences for both strands. 

Data Analyses 

COI sequences for each individual were as- 
sembled and edited using ESEE 3.0s (Cabot 
& Beckenbach, 1989). ESEE was also used 


to align our sequences with sequences ob- 
tained from GenBank and the literature. 

A distance matrix was computed using 
DNADIST of PHYLIP version 3.57 (Felsen- 
stein, 1989, 1993). Maximum likelihood trees 
were generated using DNAML (PFIYLIP 3.57). 
Twenty repetitions, with randomized input 
order and optimization by global branch re- 
arrangement, were run for each analysis. 
Bootstrap estimates (1,000 replicates) were 
made using program SEQBOOT, in conjunc- 
tion with DNAML and CONSENSE. Parsi- 
mony analyses were done using Hennig86 
(Farris, 1 988) and PAUP 3.0s with branch and 
bound searching (Swofford, 1993). In analy- 
ses using PAUP, terminal nucleotides not 
present in all sequences were trimmed from 
the analyses: using DNAML, all 645 nu- 
cleotide positions were included. 

We performed five different phylogenetic 
analyses to satisfy different conditions and 
compare among methods: (1) DNAML with 
the polyplacophoran Katharina and the rissoid 
Setia as outgroups; (2) As the polypla- 
cophoran was not much more distant from the 
ingroup taxa than was Setia, we used Setia as 
the outgroup in a maximum likelihood analy- 
sis; (3) Also with Setia as the outgroup, we ran 
Flennig86 (256 variable sites; 220 informative 
sites); (4) We subsequently added the two 
Caenogastropoda; these two taxa are better 
outgroups (i.e., within Gastropoda but outside 
Rissooidea) than either Katharina or Setia. To 
our knowledge, there are no data available for 
other rissoacean taxa that might serve as out- 
groups. For this set of taxa, we computed a 
maximum likelihood tree as follows: first, we 
empirically determined the optimal transi- 
tion/transversion ratio (ratio which minimizes 
the likelihood measure; we then used the op- 
timal ratio, 1.3 in this case, in an exhaustive 
DNAML analysis (20 iterations, global branch 
swapping); (5) we also computed parsimony 
trees for the same set of taxa, using a 1:1 
transition: transversion ratio. 


RESULTS 

Sequence alignments for 16 sequences, 
ranging in length from 578 to 645 nucleotides, 
are shown in Table 2. Individuals with identical 
sequences (414, 415; 474, 476, 477; 346, 
347; 435, 436, 439; and 453, 454) were com- 
bined for alignment and subsequent anal- 
yses. A single nucleotide at position 57 is 
missing from the Cerithium sequence; we pu- 


TABLE 2. Alignment of sequences of a 645 nucleotide fragment of the cytochrome c oxidase I gene. 
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tatively assign this deletion to a compression- 
based misreading of sequence in a GC-rich 
region. 

Felsenstein’s distance matrix of sequence 
divergence, based on data in Table 2, is given 
in Table 3. Individuals of a population differed 
by nucleotide divergence values of 0.005 ± 
0.004 (N = 5). The subspecies of Oncomela- 
nia hupensis differed by 0.148 ± 0.004 (N = 
4). The two species of Hydrobia differed by 
0.162 (range 0.161 -0.163). The genera Tric- 
ula vs. Gammatricula differed by 0.1 32 (range 
of 0.130 - 0.133). The subfamily Pomatiopsi- 
nae vs. the Triculinae differed by 0.179 ± 
0.020 (N = 12). The families Hydrobiidae vs. 
Pomatiopsidae differed by 0.267 ± 0.01 6 (N = 
21). The non-rissoacean outgroup snail taxa 
differed from the rissoaceans by 0.274 ± 
0.023. 

The unrooted maximum likelihood tree with 
the polyplacophoran as outgroup is given in 
Figure 1 . Replacement of Katharina by Setia 
as the outgroup does not change the topology 
of the remaining tree. From these trees the 
following are clear: (1) The pomatiopsid taxa 
are one distinct clade; the other rissoaceans 
form a second distinct clade; (2) The trun- 
catellids are more closely allied to the hydro- 
biids than to the pomatiopsids; (3) The rissoid 
Setia is part of the Truncatellidae-Hydrobiidae 
clade (except in Fig. 2, where it is the out- 
group); (4) The two subspecies of Oncomela' 
nia are clearly divergent; (5) The triculine taxa 
appear divergent from the pomatiopsine taxa. 
However, the triculine node is only weakly 
supported (54% bootstrapping value). 

The Hennig86 analysis using Setia as out- 
group yielded two equally parsimonious trees 
with a length of 517, a consistency index of 
0.66 and a retention index of 0.74. The Nelson 
consensus tree is shown in Figure 2. The po- 
matiopsids are in one clade, the hydrobiids 
and truncatellids in another. The results are 
the same as in the maximum likelihood analy- 
sis except that the position of Tricula is unre- 
solved: there is a trifurcation in the pomatiop- 
sid clade with Tricula not unequivocally within 
a triculine clade (it is in one of the alternative 
Hennig trees). 

The trees obtained using Stramonita and 
Cerithium as outgroups are given in Figures 3 
and 4. The maximum likelihood analysis (Fig. 
3) yields a tree for the ingroup taxa similar to 
that in Figure 1 except that Tricula is a sister 
taxon to Oncomelania rather than to Gamma- 
tricula. The PAUP-based parsimony analyses 
(Fig. 4) again clearly define two major clades. 


a Pomatiopsidae clade and a Truncatellidae- 
Hydrobiidae clade. The analysis produced 
two shortest trees of 974 steps; the set of four 
next shortest trees required one additional 
step; the latter were not included in computing 
the Adams consensus tree. The two shortest 
trees differed in their placement of Sef/'a within 
the hydrobiid clade, thus resulting in a con- 
sensus tree with a trifurcation in this clade. 


DISCUSSION 

Taxonomic decisions should not be made 
on the basis of molecular distance coefficients 
alone (Davis, 1994), but on anatomical, cyto- 
logical and developmental data within an eco- 
logical context. Thus, the COI data presented 
here must be examined in light of the avail- 
able anatomical data and the patterns of evo- 
lution evidenced in the clades shown here. 
The data provide a beginning of showing rela- 
tionships; there are insufficient hydrobiid, 
truncatellid and rissoid genera and species in- 
volved in this study to strongly support dis- 
crete rissoid, truncatellid, and hydrobiid 
clades. The data are, however, sufficient to 
answer our questions and provide the basis 
for predictions concerning family relationships 
indicated here. 

The hydrobiids and pomatiopsids are dis- 
tinct clades. The genetic data coupled with 
anatomical data reviewed in Davis (1979, 
1980, 1992) show that these clades are 
greatly divergent. The COI sequence data 
strongly support the existence of a pomatiop- 
sid clade separate from the hydrobiids. 

The distinctiveness of the pomatiopsid 
clade is further evidenced by the grouping of 
the Rissoidae and Truncatellidae with the Hy- 
drobiidae within one clade. Thus, a second 
question is answered: the Truncatellidae, by 
these data, are more allied genetically with 
the Hydrobiidae than with the Pomatiopsidae. 
However, the truncatellid branch is weakly 
supported, so additional data are required. 
Davis (1979) and Ponder (1988) hypothe- 
sized that the Truncatellidae were closely re- 
lated to the Pomatiopsidae (Rosenberg, 
1996a). The Pomatiopsidae share with the 
Truncatellidae and Assimineidae the evolu- 
tion to terrestriality from aquatic and amphibi- 
ous habitats in some clades (Rosenberg, 
1996b). However, the anatomical data avail- 
able are insufficient to resolve whether the 
Truncatellidae are closer phylogenetically to 
the Hydrobiidae or to the Pomatiopsidae. 
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Katharhm timicata 


784 



642 


Setia twricukita 474/476/477 

' Trimcatellu pulchella 480 


1000 


Trwicatella pulchella 479 

Hydrohia cf. poutieuxiui 346/347 



— Hydrohia cf. poutieuxim 351 
Hydrohia uegleda 435/436/439 


824 


— Oncomclania hupensis Impcusis 93 

1000 

Oncomclania hupensis hupensis 96 

— Oncomelania hupensis rohertsoni 45 

— Oncomelania hupensis rohertsoni 48 
(lammatricula chinensis 414 415 
Gammatricula chinensis 416 


Tricula sp. 453454 


FIG. 1. Unrooted tree based on maximum likelihood with Katharina tunicata as the outgroup. Bootstrap val- 
ues (1,000 replicates) are indicated tor each node. 


Anatomical data do indicate that the Ris- 
soidae are sister taxa to the Hydrobiidae and 
the Truncatellidae. 

A 28S ribosomal RN A phylogeny placed the 
Hydrobiidae as a sister group to the Truncatel- 
lidae and the Pomatiopsidae, with a terminal 
Tmnca/e//a- Pomatiopsidae clade divergent 
from another Truncatefla group (Rosenberg et 
al., 1997). However, as Davis (1994) pointed 
out, one must be circumspect in assessing 
phylogenies based on 28S rRNA. In the re- 
sults presented by both Davis (1994) and 
Rosenberg et al. (1 994, 1 997) involving rRNA, 
there were some considerable surprises not 
supported by comparative anatomy and al- 
lozyme-based phylogenies, due in part to the 
insufficient quantity of informative characters 
(85) to resolve a large assemblage (40) of di- 
vergent taxa. 

In the rRNA data presented by Davis (1994: 
fig. 10) there were many more differences 
among margaritiferine taxa than among all 
other unionid taxa. Cumberlandia was basal 


in the unionid clade, whereas Margaritifera 
margaritifera and M. falcatawere highly diver- 
gent from each other in a margaritiferine 
clade. The results were as if rRNA evolution 
went unexplainably berserk (greatly acceler- 
ated) in the Margaritiferidae; results not sup- 
ported by any other data (see also Ledyard et 
al., 1996). In contrast to the rRNA sequence 
data, the results with mitochondrial genes 
have been congruent with results based on 
other data, and the total weight of evidence 
leads us to consider the truncatellids within 
the same clade as, and closely related to the 
Hydrobiidae, but divergent from the Poma- 
tiopsidae. This is our hypothesis to be chal- 
lenged with results based on yet other genes. 

The COI sequence data are congruent with 
cytochrome b gene sequence data involving 
populations of Oncomelania hupensis from 
mainland China; the results showed evidence 
for two distinct subspecies: O. h. hupensis 
and O. h. rohertsoni (Spolsky et al., 1996). 
Both data sets are congruent with an al- 
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Selia tinriculata 474M76/477 

Tntncaiella piilchella 480 

Tnmcatella pulchcUa 479 

Hvdrohia nc^lccfa 435/436/439 

Hydrobiu cf. pontieitxhii 346/347 

Hydwhiti cf. pontieuxhii 351 


Tricula sp. 453/454 

Gammatricula chincnsis 4 1 4/4 1 5 

Gumwutriculu chincnsis 4 1 6 

Oncomelania hupcnsis hnpensis 96 

Oncomelania hupcnsis hupcnsis 93 

Oncomelania hupcnsis roherlsoni 45 

( )ncomclania hupcnsis robertsoni 48 


FIG. 2. Unrooted Nelson consensus tree based on Hennig86 analysis, with Setia turriculata as outgroup. 


lozyme-based phylogeny indicating the occur- 
rence of three subspecies of Oncomelania 
hupensis (the third, tangi, from Fujian Prov- 
ince, was not available for the DMA studies). 
Oncomelania h. robertsoni has a smooth 
shell, no varix, and is relatively shorter in shell 
length than O. h. hupensis, which has strong 
ribs and strong varix when living in the 
Yangtze River flood plains, but lacks ribs (al- 
though retaining a strong varix) when living 
above the effects of the annual flooding. On- 
comelania h. robertsoni lives above the Three 
Gorges of the Yangtze, whereas O. h. hupen- 
sis is found below the gorges along the 
Yangtze River drainages. The distance coeffi- 
cient between the subspecies using cy- 
tochrome b data averaged 0.110 (range of 
0.1021 - 0.1186), based on robertsoni from 
Yunnan and Sichuan and hupensis from Jiang 
Xi Province. The equivalent average distance 
in this study is 0.148 ± 0.004 (N = 4) with a 
range of 0.142 - 0.153. In this study, robert- 
soni came from Yunnan, and hupensis from 


Hubei Province, a province upstream from 
Jiang Xi.This suggests that the COI gene has 
diverged more than the cytochrome b gene 
and thus may be more useful to detect popu- 
lation and subspecies divergences and pat- 
terns of evolution. That the COI gene appears 
to have diverged more than the cytochrome b 
is of interest given the conventional under- 
standing that COI tends to be more conserva- 
tive in most taxa. 

The two species of Hydrobia differed from 
each other by a distance coefficient only about 
6% greater than the distance between the sub- 
species of Oncomelania hupensis, a surprise 
given the considerable anatomical differences 
between the hydrobiid species and the great 
geographic distance separating them (thou- 
sands of ocean miles between Bulgaria and 
Denmark). In contrast to the large differences 
between Oncomelania hupensis subspecies, 
two populations of O. h. robertsoni separated 
by over 600 km had a cytochrome b difference 
of only 0.038, a relatively small amount of dif- 
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Stnmiomta haemastoma 

( 'erithium utnUwu 

Setia furricuJafa 474/476/477 

Hydrobia cf. poiitieuxini 346/347 

Ilydrobia cf. pontieiixini 35\ 

Hydrobia ueglecta 435/436/439 

TnmcateUa pulcheUa 479 
Tnmcatclla pulchcUa 480 

— Oncomciania hupcnsis hupoisis 93 
OncowcUinia hupensi.s hupensis 96 
Oncomciania hnpcnsis robcrtsoni 48 
Oncomckmia hnpcnsis robcrtsoni 45 

Tricula sp. 453/454 

Gammatricida clhncnsis 414 415 
Gammatricnla chincnsis 4 1 6 

FIG. 3. Unrooted maximum likelihood tree using Stramonita an6 Cerithium as outgroups. Character changes 
were optimally weighted at 1.3:1 for transitionsitransversions. 


ference. It appears that the subspecies have 
diverged considerably since Oncome/an/a dis- 
persed into China about 10 million years ago 
(Davis, 1979, 1980). The considerable genetic 
divergence is also seen in allozyme data 
(Woodruff et al., 1 988; Davis et al., 1 995). The 
Nei minimum D for robertsoni vs. hupensis 
was 0.257 ± 0.077 (N = 21) (Davis et al., 
1995). Population variation within hupensis 
was 0.1 60 ± 0-085 (N = 21 ). Davis et al. (1 995) 
discuss why we do not consider these taxa as 
full species. Our question here is what has dri- 
ven such large genetic divergence when the 
snails are identical morphologically and 
anatomically (except for size and presence of 
ribs on some populations) and can replace 
each other ecologically? One hypothesis is 
that, as virtually all populations of Oncomela- 
nia hupensis throughout its range have been 
heavily infected with Schistosoma japonicum, 
coevolutionary pressures have been the 
cause. Added to this are the natural geo- 
graphical isolation of Yunnan and Sichuan 


from the central Yangtze River basin and the 
mountain ranges that separate Fujian Prov- 
ince from the other two regions. We now need 
considerable population studies of COI se- 
quences to assess the extent of population di- 
vergence across China. 

Considering the Triculinae issue, it could be 
argued that, on the basis of anatomy, the Tri- 
culinae and Pomatiopsinae are not historically 
closely related. The pomatiopsine female re- 
productive system appears considerably dif- 
ferent from that of the Triculinae, and what 
Davis has called the spermathecal duct in 
each taxon may not be homologous. The data 
presented here are congruent with allozyme 
data of Davis et al. (1994) that show general 
agreement of phylogenies based on both 
anatomical and allozyme data: Oncomelania 
hupensis is in a clade apart from the three 
triculine taxa in the allozyme study, Gamma- 
tricula chinensis, Gammatricula songi, and 
Neothcula lilii. However, the two clades are 
very close genetically. The Nei D between On- 
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Stmmonita hacmastoma 
( 'crithium atratuni 


Scliu turriculata 474/476/477 
Tnmcuielhi pulchella 480 
TrwicateUu pulchella 479 
Uyclrohia cf. pontieiixim 346/347 
Hydrohia cf. pontieuxini 35 1 
Hydrohia neglecta 435/436/439 
Oucomchmia hupensis hupensis 93 
Oncomelania hupensis hupensis 96 
( )ncomelania hupensis robertsoni 45 
Oncomelania hupensis robertsoni 48 
Tricula sp. 453/454 
Gammatricula chinensis 414/4 1 5 
Gammatricula chinensis 416 


FIG. 4. Unrooted PAUP tree: Adams consensus of two shortest trees produced by a branch and bound 
search, with transitions and transversions weighted equally. 


comelania and the three triculines was 1 .293 
± 0.412 (1.000 - 1.764); D among triculines 
was 0.890 ± 0.301 (0.689 - 1.236). Further, 
Oncomelania was less distant from Gamma- 
tricula chinensis (1.116) than was Neotricula 
lilii (1.236). Accordingly, the phenogram 
based on Nei’s D does not provide a great 
separation from the triculine taxa because of 
the closer relationship of Oncomelania to N. 
lilii. 

The Tricula of this study is a new species 
that will be described elsewhere; anatomical 
data confirm its placement in Tricula. Hen- 
nig86 produced two trees because O. h. 
robertsoni was closer to Tricula sp. (0.144) 
than it was to O. h. hupensis (0.1 50). The two 
triculine taxa differed by only 0.1 32. This very 
close relationship between Oncomelania and 
some triculine taxa shows the overall close re- 
lationship between the two pomatiopsid sub- 
families and possibly points the way to unrav- 
eling the pathway of evolution of the triculines 
from early pomatiopsids. For example, the 


Sichuan Tricula is, on the basis of anatomy, 
closely related to Yunnan species of Tricula 
and to the type species, Tricula montana from 
northern India. The average COI distance of 
Tricula sp. from O. h. hupensis is 0.188, 
whereas that of Gammatricula chinensis is 
0.195; the same set of relationships to O. h. 
robertsoni is 0.144 and 0.171. The premise 
that emerges from these data is that O. h. 
robertsoni in Yunnan, with its generalized 
morphological features, is closer genetically 
to the basal Oncomelania stock that gave rise 
to Tricula than is the derived Gammatricula 
and derived O. h. hupensis below the Three 
Gorges. This premise and initial data are con- 
sistent with the hypotheses of Davis (1979, 
1980, 1992) that pomatiopsids were intro- 
duced into Asia from the northeastern Indian 
Plate with the Himalayan orogeny, with sub- 
sequent evolution and dispersal down river 
systems. The closer genetic relationships of 
some triculine taxa to Oncomelania points the 
way to assessing the coevolution of Schisto- 
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soma spp. with both Oncomelania and Triculi- 
nae. 

All the phylogenetic analyses combined af- 
firm the monophyly of the Pomatiopsidae, al- 
though the placement of Tricula within this 
clade remains unresolved: using Katharina or 
Setia as outgroups in a maximum likelihood 
analysis, Tricula clusters with Gammatricula; 
using the closer Stramonita and Cerithium as 
outgroups, Tricula clusters with Oncomelania 
in a maximum likelihood analysis, but with 
Gammatricula in a parsimony analysis; a Nel- 
son consensus tree of Hennig86 analyses 
leaves relationships within the Pomatiopsi- 
dae as an unresolved trifurcation. In the end, 
the bootstrap value for the Tricula clade tells 
the story; placement of Tricula depends on the 
method and choice of outgroups. It is not well 
resolved within the pomatiopsid clade, but all 
evidence places it in this clade. This question 
may best be answered by increasing the num- 
ber of populations and species of Tricula and 
Gammatricula In the analyses, increasing the 
number of nucleotides analyzed, and se- 
quencing an independent nuclear gene. 

All five analyses place Truncatella in a 
clade with Hydrobia. The placement of Setia, 
however, is problematic: in analyses where it 
is not used as an outgroup, it clusters with the 
hydrobiid-truncatellid clade, but its relation- 
ships within this clade are uncertain. The un- 
resolved trifurcations are consistent with 
weak bootstrapping support for those nodes 
in the maximum likelihood analyses. 

Because the distances of KatharinaXo other 
taxa are not much different than distances 
among some ingroup members, we examined 
the number and kind of substitutions occur- 
ring in the COI gene. Substitution rates varied 
tremendously among codon positions, with 
strong constraints on amino acid changes: 
203 of the 21 5 third codon positions were vari- 
able, while 51 of 215 first codon positions and 
only ten second codon positions were vari- 
able. Translation of the nucleotide sequences 
clearly pointed out the strong constraints on 
amino acid changes in this gene: the majority 
of changes were synonymous substitutions. 
Among this group of 16 taxa, 40 amino acid 
positions were variable, of which too few (28) 
were phylogenetically informative. PAUP 
analysis of the amino acid sequences pro- 
duces 28 equally parsimonious shortest trees, 
resulting in a consensus tree with no resolu- 
tion. 

In an attempt to avoid saturation effects on 
phylogenetic tree construction, we carried out 


maximum likelihood analyses using only first 
and second codon positions: again, these 
analyses resulted in non-robust and implausi- 
ble trees with poorly supported nodes, al- 
though in this case, the pomatiopsid clade 
does hold together. We plotted pairwise tran- 
sition/transversion (Ts/Tv) ratios vs. pairwise 
distances to look for evidence of substitutional 
saturation. Although saturation must be oc- 
curring only at third codon positions, pairwise 
ratios were based on the whole sequence, as 
that is what the phylogenetic analyses were 
based on. Saturation appears to be occurring 
in pairwise comparisons with the three out- 
group taxa Katharina, Cerithium, and Stra- 
monita, with a mean Ts/Tv ratio of approxi- 
mately 1.1, and pairwise differences greater 
than about 130 nucleotides. For ingroup com- 
parisons, pairwise differences were usually 
less, varying from about 140 to 1 . Ts/Tv ratios 
were roughly correlated with nucleotide differ- 
ences, such that ratios approached saturation 
at distances greater than approximately 120 
nucleotide differences. Thus, there appears to 
be a quantitative difference between ingroup 
and outgroup taxa in the substitution level. 
The existence of a discrete ingroup provides 
further support for the robustness of the phy- 
logeny of this group. Further, although the ef- 
fect of saturation is to shorten branch length, 
in this case it is not likely to have affected the 
topology, that is, branching pattern, of the 
trees because most of the deep branches are 
sufficiently long to retain the correct branching 
pattern. 

This study should conclusively clarify two 
points: (1 ) The Hydrobiidae are not closely re- 
lated historically, biogeographically, or geneti- 
cally to the Pomatiopsidae; they are distinctly 
different families. (2) Schistosomes have 
evolved with the pomatiopsid lineage in Asia; 
hydrobiids do not transmit schistosomes. As 
implied by the term coevolved, there is a his- 
torical genetic linkage between the pomatiop- 
sids and schistosomes extending back to 
Gondwanaland. With the progression of time, 
this linkage has deepened; once lost, it can- 
not be restored (Davis, 1992). Thus, we pre- 
dict that no hydrobiid suddenly de novo, can 
become a host for any schistosome. 
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